j 

i 



i 

) 

i 

1 

NPS55-80-009 

NAVAL POSTGRADUATE SCHOOL 

! Monterey, California 




i 



i 

i 

i 

i 




AN ALGORITHM FOR POSITION DETERMINATION 
USING LORAN-C TRIPLETS WITH A 
BASIC PROGRAM FOR THE 
COMMODORE 2001 MICROCOMPUTER 

by 

R. H. Shudde 

/) 

March 1980 



Approved for public release; distribution unlimited. 
Prepared for: 



S 



i 

i 

f 



Naval Research 

FEDDOCS , Virginia 22217 

D 208.1 4/2:NPS-55-80-009 




NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 



CUDl£ v KNOX LIBRARY 
imi POSTGRADUATE SCHOU 

MnHTtREY. CA 93940 



Rear Admiral J. J. Ekelund 
Superintendent 



Jack R. Borsting 
Provost 



The work reported herein was supported by the Office of Naval Research, 
Fleet Acitivity Support Division, Code 230 and the Commander, Patrol Wings 
Pacific as part of the Tactical Development and Evaluation Program. 

Reproduction of all or part of this report is authorized. 



UNCLASSIFIED 



.!i)LL N rJNUA uonwii 
»VAl TO ST GRAD I ATE cr 



SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. report number 
NPS55-80-009 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE (and Subtitle) 

AN ALGORITHM FOR POSITION DETERMINATION USING 
LORAN-C TRIPLETS WITH A BASIC PROGRAM FOR THE 
COMMODORE 2001 MICROMPUTER 


S. TYPE OF REPORT & PERIOD COVERED 

Technical 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHORO) 

R. H. Shudde 


0. CONTRACT OR GRANT N UM 8 E R(e) 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, Ca. 93940 


10. program element, project, task 

AREA S WORK UNIT NUMBERS 

65155N, R0130 ,NR 128-217 
N0001479WR90234 


11. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, CA 93940 


12. REPORT DATE 

March 1980 


13. NUMBER OF PAGES 


14 MONITORING AGENCY NAME 4 ADDRESS^// different from Controlling Otttce) 


15 SECURITY CLASS, (of thle report) 

Unclassified 


15a. DECLASSI FI CATION/ DOWN GRADING 
SCHEDULE 


16. DISTRIBUTION STATEMENT (ot this Report) 

Approved for public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (of the abetract entered In Block 20, it different from Report) 


10. supplementary notes 


19. KEY WORDS (Continue on reverse eide if necessary and identify by block number) 

Lor an. Fixing, Geodetic Distances 

Loran-C, Hyperbolic Fixing, 

Navigation, Radio Positioning, 

Position determination, Geodetics, 


20. ABSTRACT (Continue on reverae eide if necessary and identify by block number) 

An algorithm for position determination using Loran-C triplets and an imple- 
mentation of the algorithm in the BASIC language of the Commodore 2001 micro- 
computer are presented. Modifications of the geodesic direct solution and 
reverse solution algorithms of P. D. Thomas are also presented; these modi- 
fications eliminate all of the special positioning and quadrant determina- 
tion requirements of the original algorithms. 



DD 1 jan^73 1473 EDITION OF 1 NOV 65 IS OBSOLETE 

S/N 0 102-014- 6601 i UNCLASSIFIED 

SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 



AN ALGORITHM FOR POSITION DETERMINATION 



USING LORAN-C TRIPLETS WITH A 
BASIC PROGRAM FOR THE 
COMMODORE 2001 MICROCOMPUTER 



by 

R* H. Shudde 

Naval Postgraduate School 
Monterey, CA 93940 



March 1980 



The program in this report 
is for use within the Navy 
and is presented without 
representation or warranty 
of any kind. 



CONTENTS 



Pace 



A. Introduction 1 

B. Loran-C Fixing Algorithms 2 

C. The Direct Solution Algorithm 10 

D. The Reverse (Inverse) Solution Algorithm 12 

E. Program Accuracy 14 

F. References 17 

Appendix: BASIC Program Listing 18 



ABSTRACT 



An algorithm for position determination using Loran-C trip- 
lets and an implementation of the algorithm in the BASIC 
language of the Commodore 2001 microcomputer are presented. 
Modifications of the geodesic direct solution and reverse 
solution algorithms of P. D. Thomas are also presented; 
these modifications eliminate all of the special position- 
ing and quadrant determination requirements of the original 
algorithms . 



A. 



Introduction 



Loran-C position determination programs for the Hewlett- 
Packard HP-67 programmable calculator for use aboard U. S. 

Navy patrol aircraft (P-3) have been prepared and are reported 
elsewhere [Ref. 1]. During the development and testing of 
the HP-67 programs the need for a high accuracy program for 
cross checking purposes arose. The program presented here was 
developed as a refinement of the HP-67 algorithm. 

The basic methodology of this program is contained in 
the formulas of Paul D. Thomas [Ref. 2] for the direct and 
reverse problems of oblate spheroidal geodesy. The direct 
problem is to determine the latitude and longitude of a point 
P 2 when the latitude and longitude of a point P-^ as well as 
the azimuth and distance of P 2 from P^ are known. The 
reverse (or inverse) problem is to determine the distance, 
forward azimuth, and backward azimuth between two points, P-^ 
and P 2 , when the latitude and longitude of both P^ and P 2 
are known. In both problems P-^ and P 2 are on the surface of 
a spheroid. Unfortunately, Thomas' algorithms contain several 
restrictions and numerous subcases that are too awkward to be 
programmed on the HP-67. For this reason portions of these 
algorithms were rewritten so as to eliminate all of the original 
restrictions and subcases. These modified algorithms for the 
direct and reverse problems are presented here with notations 
about the modifications made. 
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B . Loran-C Fixing Algorithms 

The principles of Loran lines of position (LOP's) and 
fixing are adequately covered in Reference 3 and will not be 
repeated here. 

The basic Loran-C equation [Ref. 4] can be written as 

T = [Tg + p (Tg) ] - [T m + p(T M ) ] + [T fi + p (Tg) ] + 6 (1) 

where 

T is the "indicated time difference" in microseconds, 

T m is the distance, in microseconds, from the master 
to the receiver, 

Tg is the distance, in microseconds, from the slave 
to the receiver, 

Tg is the distance, in microseconds, between the 
master and the slave, 

6 is the assigned coding delay, in microseconds, 
and p(T) is the secondary phase correction, in micro- 
seconds, for a surface water path of length T. 

The quantity 

At = [ Tg+p (Tg) ] + 6 

is a constant for each master/slave pair. The following 
World Geodetic System 1972 (WGS 72) values have been adopted 
for Loran-C navigation [Ref. 4] : 

Vq = 299792458 meters/second is the velocity of 
light in free space, 

ri = 1.000338 is the index of refraction of the 

surface of the earth for standard atmosphere 
and 100 kHz electromagnetic waves. 
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a e = 6378135.000 meters is the equatorial radius 
of the earth 

and f = 1/298.26 is the flattening factor (l-b/a e , 

where b is the polar radius) of the earth. 

With these parameters the secondary phase correction for 
an all seawater path has been taken to be of the form 

p(T) = a Q /T + a 1 + a 2 T (2) 

where T is in microseconds and, 

1. For T>537 ysec: 

a Q = 129.04398, 
a ± = -0.40758, 
and a 2 = 0.00064576438. 

2. For t< 537 ysec: 

a Q = 2.7412979, 
a ± = -0.011402, 
and a 2 = 0.00032774624. 

If one uses a spherical approximation to the earth's 
surface (f = 0) , then a spherical hyperbola can be represented 
by the equation [Ref. 3, page 175] 

^ g pj £ _ 2 cos 2 a — cos 2c ( 3 ) 

sin 2c cos co + £ sin 2a 

where the origin of the coordinate system is at the prime 
focus of the spherical hyperbola, 2c is the spherical arc 
joining the foci, 2a is a constant for any one hyperbola, 
and r and co are the spherical coordinates of a point on the 
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hyperbola. If the base line of the coordinate system is the 
arc joining the foci then w is the spherical polar angle from 
the baseline to a point P on the spherical hyperbola and r is 
spherical polar distance (or arc) from the prime focus to P. 

Using the Loran system we take Z, - +1 if the prime focus is 
at a master station and we take C = -1 if the prime focus is 
at a slave station. 

If we take v = Vq/ti to the velocity of 100kHz 
electromagnetic radiation at the earth's surface then, for a 
spherical earth, we can relate the parameters in Equations 1 
and 3 as follows: 

2c - vr B /a e , 

and 2a = v(Tg - T H )/a e - 

Using the spherical approximation for now, we see that 2c is 
known for any Loran pair. The "indicated time delay" T is 
measured by the receiver at point P, and to determine a hyper- 
bolic line of position we must determine 2a, but T c - t m 
cannot be computed from Equations 1 and 2. If a^ were zero 

in Equation 2, then it would be possible to determine T„ - T„ 

b M 

uniquely. As a first approximation we use the following 
parameters in Equation 2: 

y a 0 " °' 

a-L = -0.321, 

and a. 2 = 0.000635. 

These values have been obtained by setting ag = 0 and determining 
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a-^ and a 2 by linear regression of the T > 537 values over 
the interval of 1000 < T < 8000. This approximation is quite 
good (within 0.03 ys) for distances up to 10,000 microseconds 
where small changes in the LOP's can cause large position 
errors. At short distances the error increases from 0.05 ys at 
1000 ys to 0.58 ys at 10 ys; although these errors are large for 
small distances, the LOP's are not as sensitive to these 
changes as they would be at large distances. These errors are 
well within the 4 n.mi (24 ys) accuracy that was requested by 
COMPATWINGSPAC . When this approximation is substituted into 
Equation 1, we obtain 



[T s + + a 2 T s ] - [t m + = T - At, 



or 



T s " t m = (T - At)/ (1 + a 2 ) 



(4) 



and hence 2a = v(T c -T M )/a is determined for use in the 

s M e 

spherical approximation. 

Consider a Loran-C triplet with the master stations colocated. 
Let and £ 2 denote the azimuth angles of slave 1 (S^) and 
slave 2 (S 2 ) , respectively, measured from North toward the East 
from the master stations (M) (see Figure 1). Further, let a 
and r be the azimuth and spherical polar arc (distance) of the 
receiver (R) from M. For this geometry. Equation 3 can be 
written in the form 



tan r ^ 



B 



COS (a - ^i) + A i 



(5) 



5 



N 




Figure 1. Geometry of a Loran Triplet and a Receiver. 
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where 



A i = sin 2a., 



= cos 2a^ - cos 2c^ , 
and = sin 2c^ 

for the i th Loran pair, i = 1,2. Since r-^ = r 2 = r, 
we can eliminate tan r between the two equations. The 
resulting equation can be rewritten as 



C cos a + S sin a = K, (6) 

where 





C " B 1 C 2 


cos ^2 ~ 


B 2 C 1 


COS 


5 i< 




S = BjC 2 


sin K 2 - 


B 2 C 1 


sin 


r i' 


and 


K - B 2 A 1 


b 1 A 2* 









If we define p > 0 and y by the equations 

p cos y = C , (7) 

and p sin y = S , 

then 



p = /c 2 + S 2 , 

and y = qatn (S,C) . 

Here the function qatn(y,x) is the arctangent of y/x adjusted 
for the proper quadrant according to the signs of x and y. 
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A compact form of this function is 



_1 * 

qatn(y,x) = tan ^ + Tft(x < 0?) 

x + 10 t (x = 0?) 

where t(z) =1 when z is true 

and t ( z ) =0 when z is false. 



When convenient we will use the notation qatn(y/x) inter- 
changeably with qatn(y,x). Now we can substitute Eq.(7) into 
Eq. (6) and solve for 

a = y ± cos“ 1 (K/p) (8) 



to obtain the azimuth angle a of the two points of inter- 
section of the spherical hyperbolic LOP's. Finally we can 
obtain a value for r by substituting each a into either 

Eq. (5) . We find that 



r = qatn 



*i 

C ^ cos (a-£ ^ ) + A i 



for i = 1 or 2, 



To obtain the first solution, r and a are entered into 
the "direct" solution algorithm; latitude cf> and longitude A 
are the outputs. 

The solution is improved using a two-variable Newton- 
Raphson search for the zeroes of the functions 



f i ' ^) = f i = ( d si - d M )/ a e " 2a i for 1 = 1 ' 2 ' 

where d M , d g ^, and d g2 are the distances from R to M, 
S^, and S 2 , respectively, and are computed using the 
"inverse" solution algorithm. The flow of the improvement 
algorithm follows: 
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1. Input the latest approximation to the receiver latitude 
and longitude A. 

2. Compute d^, d g ^ and d g 2 using the "inverse" solution 
algorithm . 

3. Compute T M = d M /v and T g ^ = dg^/v for i = 1,2. 

4. Compute T gi - T M = (T-At). - p(T gi ) + p(T M ) and 

2a i = v (T si" T M )/ a e for i = 1,2m 

5. Compute f.(<f>,A) = (d c .-d M )/a - 2a- for i = 1,2. 

j. o x w 6 l 

6. Compute 

i _ f A ( 4» + A*, A) - f ± ( 4>, X) 

3(f) A c)> 

and 

iii _ +&\ i - for i . 1<2 . 

3 A ax 



- 4 

(AA = Acf) = 10 radians is used in the BASIC program.) 
7. Replace <J> by 




and replace A by 

A 

where 

j = _!£i ^2 _ 

3 4> 3 A 3 A 34> 

8. Repeat from Step 1 until <f> and \ are stationary. 
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c. 



The Direct Solution Algorithm 



This direct solution algorithm is a modification of the 
second order in flattening (f) algorithm given by Thomas 
[Ref. 2, pp. 7-8]. Thomas' notation has been followed as 
closely as possible for ease of comparison of the algorithms. 
The qatn function is defined in the previous section. West 
longitudes and South latitudes are negative. We are given the 
point P 1 (c() 1 ,A 1 ) on the spheroid, where <j> lf are the geodetic 
latitude and longitude (geographic coordinates); the forward 



from these we are to find the geographic coordinates ‘t 1 2 ' ^ 2 and 
the back azimuth c^-^. The given quantities are <j>^, X a-^ 
and S. No assumptions about the relative location of and P 
are required. The modified direct solution algorithm is: 




0^ = tan 1 [(l-f) tan cj>^], M 



cos 0^ sin ct^ 2 ' 



N = cos 0^ cos ct^ 2 » c \ 



fM, c 2 = f ( 1-M 2 ) / 4 



D = (1 - c 2 ) (1 - c 2 - CjM) , P = c 2 ll + (l/2) Cl M]/D 



= qatn(N, sin 0 ^), d = S/(a e D) 



u = 2(a.-d), W = 1-2P cos u, V = cos (u+d) 



X = c 2 2 sin d cos d (2V 2 -1) , Y = 2PVW sin d 



A a =d+X-Y, Za= 2a - Aa 



“21 



qatn [-M, - (N cos Ao - sin 0^ sin Ao) ] , 



K = (1-f) [M 2 + (N cos A 0 



a >2,1/2 

- sm 0^ sm Act) ] , 

<J >2 = tan '*'[(sin 0^ cos Ao + N sin Aa)/K] , 

An = qatn(sin Ao sin a^( c os 0^ cos Aa - sin 0^ sin Ao cos a^) • 
H = c^(1-C2)Ao - c^C 2 sin Ao cos La, 

AA = An - H, A 2 = A 1 + AX 

In addition to the introduction of the qatn function for 
proper quadrant determination the following changes have been 
made to the original algorithm: 

1. 9g is no longer computed from the equation M = cos 0g. 
Since P 1 is no longer required to be westerly of P 2 , 

^12 * s no l° n 9 er required to lie in the interval 

[0, 180°]. Consequently M, which was originally 
required to be positive, can be negative. With this 
change Gg was no longer properly determined from 
M = cos Gq. 

2. With the elimination of Gg the determination of 
o 1 f rom the equation cos = cos 0^/sin Gg became 
impossible. The following figure shows the spherical 
triangle involving and Gg. 

9O°-0g 

1 
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Using the "4-parts formula" [Ref. 5, pg . 12] we can 
determine directly from the equation 

= qatn(cos 0^ cos a^ 2 , s i n ®]_) 

3. The original equation for <|> 2 r 

tan 4 > 2 = -(sin 6^ cos Ao + n sin Ao) sin a 2 ^/(l-f)M 

becomes indeterminate when a^ 2 = 0° or 180° since 
both M and sin become zero simultaneously. 

From the equation for a 21 one can determine that 

sin a 21 = -M/[M 2 + (N cos Ao - sin 0^ sin Aa) 2 ]-*-/ 2 . 

Using this equation to eliminate sin a 2 in the equation 
for tan <f> 2 allows the M in the numerator and denominator 
to cancel thus eliminating the indeterminancy in the 
equation for <l> 2 . 

With these changes the restrictions that be west of P 2 , 
and that if P 1 and P, both have negative latitude the symmetric 
positive latitude problem be solved, have been eliminated. 

Also eliminated are rules for the determination of the quadrants 
of a 2 ^ and An. 

D . The Reverse (Inverse) Solution Algorithm 

This reverse solution algorithm is a modification of the 
second order in flattening (f) algorithm given by Thomas [Ref. 2, 
pp. 8-10]. Thomas' notation has been followed as closely as 
possible for ease of comparison of the algorithms. The qatn 
function is defined in Section B. West longitudes (X) and 
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South latitudes ( <j>) are negative. We are given the points 
P 2^2'^2^ on t ^ ie s P^ lero: *- < 3 and are to find the 
distance S between the points and the forward and back azimuths, 
0^2 and o^. Given quantities are <f>^ , A-^, 4> 2 and X 2 * No 
assumptions about the relative location of and P 2 are 
required. The modified reverse solution algorithm is: 

= tan -1 t(l-f) tan 4 k], i = 1,2, 

9 m = ( 0 i +e 2 )/ 2 ' A0 m = ( 0 2 - 0 i)/ 2 ' AX = X 2 - X x , 

2 2 2 2 

AX = AX/2, H = cos A0 -sin 0 - cos 0 - sin A0 = cos 0. cos 0_ 

m m m m m 12 

L = sin 2 A0 m + H sin 2 AA m = sin 2 (d/2), 

1-L = cos 2 (d/2) , d = cos _1 (l-2L), 

U = 2 sin 2 0 m cos 2 A0 m /(1-L), V = 2 sin 2 A0 m cos 2 0 m /L, 

X=U+V, Y = U - V, T = d/sin d, D = 4T 2 , 

E = 2 cos d, A = DE, B = 2D , C = T - (A-E)/2 
n x = X ( A+CX ) , n 2 = Y (B+EY ) , n 3 = DXY , 

6 ]d = f (TX-Y ) /4 , 5 2 d = f 2 (n 1 ~n 2 +n 3 ) /64 , 

S = a e (T- 6 1 d+6 2 d) sin d, F = 2Y - E(4-x) 

M = 32T - ( 20T-A) X - (B+4)Y, 

G = fT/2 + f 2 M/ 6 4 , Q = -(FG tan AA)/4, 
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A V = ( AX + Q)/ 2 ' 

= qatn (-sin AO^ cos A A^, cos ^ sin A^ ' ) » 
t 2 = qatn(cos A6 m cos A^', sin ^ sin A^')r 

°12 = t l + t 2' a 21 = fc l “ t 2 



The only changes made to the original algorithm are the 
computation of t-, and t 2 and the determination of and 

as the sum and difference of t-^ and t 2 » Minor though these 
changes may seem, they have eliminated the requirement that 
P^ is West of P 2 and they have eliminated four cases for the 
quadrant determination of and a 21‘ Quadrant determination 

is still required, but it is implemented with the qatn 
function; the qatn function is available on the better handheld 
calculators in the form of the rectangular-to-polar function. 



E . Program Accuracy 

Several programs were written for the Commodore 2001 
computer. One program implemented only the direct and reverse 
solution algorithms presented in the previous two sections. 
Initially these algorithms were validated against the long 
line computations with the Clarke 1866 spheroid model presented 
in Reference 2. Later the reverse solution algorithm was 
compared with the WGS 72 data for 40 Loran-C pairs contained in 
Reference 4 (data for the 9930 stations were not included) . 

In these comparisons the algorithm error for the baseline 
distances was 0.16 meters high for the 7980W pair and 0.15 
meters low for the 7980Y pair. For the remaining stations 
all of the baseline errors were high, but the largest error 
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was only 0.08 meters. It should be mentioned that the Micro- 
soft BASIC for the Commodore 2001 carries 10 digits of floating 
point internally, but will display only a maximum of nine 
digits; the baselines reported in Reference 4 are printed to 
the nearest 0.0001 meters. The computed one way baseline time 
plus secondary phase correction was also computed and compared 
to the 40 Loran-C pairs data. in three cases the algorithm 
was high by 0.01 ys, in three cases the algorithm was low by 
0.01 ys , and there was no error in the remaining cases. 

Another program was used to generate the indicated time 
delay between Loran-C pairs and test positions for a receiver. 
Combining the time delay for Loran-C pairs into allowable 
triplets enabled the position fixing program in the Appendix 
to be evaluated. Typical results are presented in Table I. 

In all cases the Commodore 2001, running at BASIC interpreter 
speed required 8 seconds from entry of the time delay input to 
the display of the first pair of solutions. The improvement 
algorithm will improve only one user selected solution. Each 
iteration requires 14.4 seconds, and each improvement is dis- 
played so that the user can monitor the improvement process. 
Usually the solution has stabilized after two iterations, but 
occasionally a third iteration is generated. in Table I the worst 
fix took six iterations for the triplet comprised of pairs 9930W 
and 9930X where the radial position error is 0.91 n.mi. In fair- 
ness to the algorithm it should be noted that the three stations 
are 776 n.mi, 2130 n.mi, and 2514 n.mi from the receiver and 
are all within 6° of azimuth as viewed from the receiver. 
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Table I 



Receiver Location: 52°N, 35°W 



Station 

Pair 


Computed 

Delay 


Station 

Pair 


Computed 

Delay 


7970W 


32845.12 


7930W 


16673.77 


7970X 


18889.03 


7930X 


30899.81 


9930W 


16067.09 


7930Z 


49369.30 


9930X 


28012.61 






Triad 


First 

Solution 


Improved 

Solution 



7970W/7970X 


52°00 ' 
34°59 ' 


14 "N 
34"W 


9930W/9930X 


51°56 ' 
35 21' 


13"N 
23 "W 


7930W/7930X 


52°00 ' 
35°00 ' 


01"N 

01"W 


7930W/7930Z 


52°00 ' 
35°00 1 


01"N 

01"W 


7930X/7930Z 


52°00' 
35 00' 


01"N 
01 "W 


7930X/7970W 


52°00 ' 
34°59 ' 


04"N 

57"W 


7930Z/9930X 


51°59' 
34 59' 


52"N 

26"W 



52°00'00".14N, 

34°59 1 59" . 66W 2 

51°59 ' 45" . 25N , 

35°01 ' 25" . 37w 6 

51°59 1 59" . 69N , 

35°00 1 00" . 02W 2 

5 2° 0 0 ' 0 0 " . 0 3N , 

35°00 ' 00" . 07W 2 

52°00 ' 00" . 00N , 

34°59 ' 59" . 96W 2 

52°00 ' 00" . 01N , 

34°59 ' 59" . 96W 2 

51°59 ' 59" . 60N 

34°59 ' 58" . 27W 2 
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APPENDIX. BASIC Program Listing 

This appendix contains the listing of the BASIC position 
fixing program for the Commodore 2001. A number of the 
cursor control and formatting symbols on the listing require 
explanation. These symbols are given below by line number at 
first occurrence. 



130 


Reverse 


f ield 


'heart' means 'clear the screen' 


130 


Reverse 


f ield 


'Q' 


means 


'move cursor down' 


250 


Reverse 


f ield 


'R' 


means 


'turn on field reverse' 


250 


Reverse 


field 


1 1 


means 


'turn off field reverse' 


250 


Use of 


the field 


reverse 


symbols highlights the 'Y' 




'YES' and the 


' N ' 


in 'NO 


' to prompt the user to type 




'Y' or 


'N' on 


the 


keyboard. Similar highlighting is 




used on 


lines 


390 


, 790, 


800 and 810 



740 CHR$(34) generates the double quote symbol to designate 

'seconds of arc' 

740 The symbol between the double quotes which are between 
N$ ( 1) and N$(2) is used to designate '°', the 'degrees 
of arc' symbol 

910 The first symbol following the double quote means 'move 
cursor up', and the next symbols each mean 'move cursor 
right ' 

The DATA statements on lines 50001 through 50044 contain 
the pertinent parameters for each Loran-C pair. The fields 
contain : 
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1 . 



The Loran-C station pair designator. 

2. At, the sum of the coding delay plus one way base - 
line time, including the secondary phase correction 
for an all seawater path, in microseconds. 

3. The master station latitude. 

4. The master station longitude. 

5. The slave station latitude. 

6. The slave station longitude. 

Negative longitudes are West longitudes. The latitudes and 
longitudes appear to be in decimal form, but the actual format 
is DDD.MMSSFF where 



DDD 


designates 


degrees , 


MM 


designates 


minutes , 


SS 


designates 


seconds , 


FF 


designates 


hundredths of seconds 



19 



50001 

50002 

50003 

50004 

50005 
5000b 
50007 
50003 
50003 
50010 
5001 1 

50012 

50013 

50014 

50015 
500 1 6 
50017 
50013 
50013 
50020 

10021 
10022 
.10023 

50024 

50025 

50026 

50027 
50023 
53023 
50030 
5303 1 

50032 

50033 
50334 
50035 
5003b 
3O037 
53033 
50033 

50040 

50041 
53042 
50G43 

59044 



cr 

5 

5 



IiftTfl 

HATH 

hath 

DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DAT A 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DhT a 
DAT A 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 



4390V 


15972 


> cLZ* . 


1 6 . 


444395, 


-163. 


303120, 


' 20. 


144316, 


-155. 


538970 


4330V . 


34253 


. 18. 


16. 


444395, 


— 1 69. 


303 1 20 ■ 


ZO « 


234177, 


-178. 


1 73820 


5330:4. 


13131 


O O 

a 1 L 1 . 


46. 


4 


82720, 


-067 . 


553771, 


41. 


151193, 


-063. 


583303 


5330V. 


28755 


.02, 


46. 


4 


82720, 


—067 . 


553771, 


46. 


463218, 


-053. 


102816 


5330:4 . 


13343 


. 60. 


51. 


5 


75878, 


-122. 


220224, 


55. 


262085, 


-131. 


151365 


5330V. 


28927 


. 36 . 


51. 


5 


75878, 


-122. 


220224, 


■ 47. 


0-34799 , 


-113. 


443953 


5330Z. 


42266 


. 63 . 


51. 


cr 


-7 err . -70 

1 • J 1 O ; 


-122. 


220224 , 


. 50 . 


362972, 


-127. 


212935 


7330 W. 


15068 


. 02. 


59. 


5 


9 1 727 , 


-045. 


102747, 


■ 64 . 


542658, 


-023. 


552175 


733054. 


27803 


• r' i - ' . 


59. 


5 


91727, 


-045. 


102747, 


' 62 . 


175968, 


-007 . 


042671 


7330Z. 


43212 


. 20 . 


53. 


cr 


31727, 


-045. 


102747, 


1 46. 


463218, 


-853. 


102816 


7360:4 . 


13804 


.45. 


63 . 


i 


34281, 


-142. 


483 1 90 , 


cr-? 

1 Jl a 


26202 1 , 


-152. 


221122 


7360V . 


29651 


. 14. 


63 a 


1 


34281 , 


-142. 


483130, 


.55. 


262 085 j 


-131. 


151965 


7370 W. 


30965 


.64. 


6iz! a 


1 


75368, 


-007. 


042671 , 


• 54 a 


432980 , 


, +003. 


173633 


7379:4 . 


15043 


. 10. 


62 a 


1 


75368, 


—00 r . 


042671 , 


• 68 a 


3806 1 5 , 


* +0 1 4 . 


274780 


7370V. 


48944 


. 53. 


62. 


X 


75368, 


-007. 


04367 1 


.64. 


542658 , 


. -023. 


552175 


7370Z. 


632 1 6 


. 30. 


u. a 


1 


75968, 


-007 . 


0426 r 1 . 


.70. 


545261, 


. "008. 


435863 


7380 W . 


12803 


.54. 


39. 


5 


33874, 


-035. 


1 00330 . 


• 30 a 


433382, 


■ -030. 


434368 


7380K . 


2!7 44 3 


"1 i”i 


30. 


5 


33074 , 


-085. 


1 00330 . 


• 26. 


315581 , 


. -097. 


500003 


7380V. 


45201 


a C*0 


30. 


5 


33874, 


-085 . 


1 00330 . 


1 

'dim 


015849, 


• -080 . 


065352 


7980Z 


.61542 


> l‘ ul J 


. 36. 


5 


93374, 


-085. 


1 00930 , 


. 34 . 


034604 , 


.-077. 


54467 6 


7330:4 .. 


12755 


.97; 


» 38 . 


5 


22061 , 


016. 


430596 . 


. 35 . 


3 l 2888 . 


. 012. 


3 1 2336 


7390V.. 


y ci. d r* 0 


. 30 ■ 


0 0 


5 


22061 , 


016. 


438596 , 


. 48. 


582035 , 


. 027. 


528152 


7390Z.. 


50993 


. 69; 


1 yy . 


cr 

•J 


22061 , 


01 6 . 


430596 , 


' 42. 


033643, 


. 003. 


121538 


3978W, 


i 4355 


all; 


. 39. 


5 


10754, 


— 08r* . 


291214, 


. 30. 


533874 , 


—085. 


1 00330 


b'37o:4.. 


31 162 


. 06; 


39 a 


cr 

j 


10754, 


-087 . 


291214, 


> 42 . 


425068 , 


. -876. 


493386 


8970V, 


47753 


a 7 4 ; 


' 39 a 


cr 


1 0754 , 


-087. 


231214. 


.48. 


364 3 y 4 , 


. -094. 


331847 


9930k! , 


13695 


.51, 


34 ! 


0 


34634 , 


“* y f * r* a 


54467 6 . 


' Cl l m 


015849.. 


, -0O0 . 


065352 


3330:4 , 


36389 


a 66, 


34. 


fi 


34604 , 


-077. 


54467 6 , 


> 46. 


463218, 


. —053 . 


102816 


9330V , 


5254 1 


a 31 , 


1 34. 


0 


34604 , 


— 0 1 '' (' • 


54467 6 , 


.41. 


151133, 


-863. 


583909 


3930 Z , 


68560 


— !* •“ 1 
. 1 C- , 


34. 


0 


34604, 


-077. 


544676, 


. '“ ,C 1 


510754, 


.-037. 


291214 


934y W , 


1 3 7 3 6 


a 90, 


, 39 . 




30662 , 


-118. 


495637 , 


.47. 


034793, 


.-113. 


443353 


3940:4 , 


28834 


a 50; 


■ 33 . 


3 


30662 , 


-118. 


495637, 


1 -Z' C' ■ 


465633 , 


.-122. 


294453 


3940V , 


4 1 367 


a 30; 


, 39 . 




30662 , 


-118. 


495637 , 


-icr 
• •-.» > J • 


191318, 


.-114. 


481743 


9360 IT, 


1 3737 


a 20 ; 


42. 


4 


25060 , 


— 0 1' 6 . 


493.386 , 


.46. 


482728 , 


. -067. 


553771 


3360:4, 


26363 


a 93; 


42. 


4 


25060, 


-076. 


4 9 3 3 8 6 , 


. 41 . 


151193, 


- -063. 


583983 
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42221 


a 65; 


42. 
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-076. 


493386, 


. 34 . 


034684 , 


. — 0 r r . 


544676 


3360Z , 


57162 


a 06, 


42. 


4 


25060, 


-076 . 


433386 , 


. 33 . 


518754, 


. -087. 


291214 


9370 N , 


15283 


a 94 ; 


24. 


4 


8041 , 


141. 


13230 , 


• 24. 


17877 , 


. 153. 


58515 


9370:4 , 


36685 


4 ”1 

a i 2- ; 


1 24. 


4 


804 1 , 


141. 


1 3290 . 


■ 42 . 


443700, 


. 143. 


430986 


z^ y »• 1.1 t , 


59463 


^ .”1 

a .i ; 


24. 


4 


8041 , 


141. 


13238 . 


■ 2'6 . 


36*2439 , 


1 1 . 


08562 1 


397GZ, 


80746 


- 79 , 


24. 


4 


8041 , 


141. 


19290 


.09. 


324566. 


. 1 38 . 


095523 


999 LI 1 "! , 


1 4875 


a 1 u. J 


cr ~ 1 

l' a 


0 


30988, 


-170. 


145381 


.52. 


494585. 


.+173. 


185231 


9990V , 


32063 


. 09 , 


' 57. 


0 


30388 , 


-170. 


145981, 


■ 65. 


144812, 


.-166. 


531447 


993 0Z , 


46530 


. 10 , 


cr~? 

' •-« . a 


0 


90988 , 


-170. 


145981 


.57. 


262021 , 


.-152. 


221122 
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Coverage of Loran-C Systems 
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Location 
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Central Pacific 
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East Coast, Canada 


5990 


West Coast,' Canada 


7930 


North Atlantic 


7960 


Gulf of Alaska 


7970 


Norwegian Sea 


7980 


Southeast U.S.A. 


7990 


Mediterranean Sea 


8970 


Great Lakes 


9930 


East Coast, U.S.A. 


9940 


West Coast, U.S.A. 


9960 


Northeast U.S.A. 


9970 


Northwest Pacific 


9990 


North Pacific 
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